Fast cubature of volume potentials over rectangular domains 
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Abstract 

In the present paper we study high-order cubature formulas for the computation of 
advection-diffusion potentials over boxes. By using the basis functions introduced in the 
theory of approximate approximations, the cubature of a potential is reduced to the quadra- 
ture of one dimensional integrals. For densities with separated approximation, we derive a 
tensor product representation of the integral operator which admits efficient cubature pro- 
cedures in very high dimensions. Numerical tests show that these formulas are accurate and 
provide approximation of order 0(h 6 ) up to dimension 10 s . 

Keywords. Multi-dimensional convolution; Advection-diffusion potential; Tensor product 
representation; Higher dimensions 
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1 Introduction 



High-dimensional volume potentials arise in many mathematical models in the field of physics, 
chemistry, biology, financial mathematics and many others. In recent years, tensor product 
approximation has been recognized as a successful tool to overcome the "curse of dimensionality" 
and treat high-dimensional integral operators as described, for example, in [3, 4, 6, 2]. 

In the present paper we propose to combine high-order semi-analytic cubature formulas, 
obtained by using the method of approximate approximations (see [11] and the reference therein), 
with tensor product approximations. 

Cubature formulas based on approximate approximations for volume potentials over M n and 
over bounded domains have been considered in [10] and [9], respectively (see also [11])- The 
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cubature of high-dimensional volume potentials over the full space and over half-spaces has been 
studied in [7] and [8]. Now we consider the volume potential 



£a/(x)=/ K A (x-y)/(y)dy, 
J[P,Q] 



(1.1) 



with the fundamental solution 



«a(x) 



(2tt)"/ 2 V A 



xKl-n/2 



^„/ 2 _i(A|x|),AGC\(-oo,0], 



over rectangular domains [P,Q] = YYj=i[PjiQj] c Here -ftTj, is the modified Bessel function 
of the second kind (see [1, 9.6, p. 374]). 

The function u = Kf provides a solution of the modified Helmholtz equation 



(-A + X 2 )u= { 



{ /(x) XG [P,Q] 
otherwise. 



For A = 0, then 



' 1 i 1 

log 



K (x) = < 



2vr ' |x| 

r(f-i) i 

4 ?r n/2 



X 



n-2 



n = 2, 
, n > 3 



is the fundamental solution of the Laplacian. 

The theory of approximate approximations proposes semi-analytic cubature formulas for 
volume potentials by using quasi-interpolation of the density / by functions for which the integral 
operator can be taken analytically. Approximate quasi-interpolant has the form 



x-/w^-" /2 E/(^(^) 



where h and V are positive parameters and rj is a smooth and rapidly decaying function which 
satisfies the moment conditions of order N 



/ t/(x) x°dx = <? 0>a , < |a| < N. 
If / € C^(R n ), it is known ([11]) that 

|/(x) - M h ,vf(x)\ < c(Vvh) N \\V N f\\ L °° + e k (Vvh) k \V k f( x )\ 

with 



(1.2) 



k=0 



e k < V \Vk^r](VVm)\; lim V |V fc ^(VPm)| = 0. 

meZ"\{0} meZ™\{0} 
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If we replace / in (1.1) by the quasi-interpolant 

V-/* £ f(hm)J^) (1.3) 

we don't obtain good approximations because (1.3) approximates / only in a subdomain of 
[P, Q] with positive distance from the boundary. To avoid this difficulty we extend / with 
preserved smoothness in a larger domain. Obviously the quasi-interpolant of the continuation 
/ approximates / in [P, Q] . Assume that there exists C > such that 

ll/llw^(R") - C Wf\\w£([P,Q])- 

Since rj is a smooth and rapidly decaying function, for any error e > one can fix r > and the 
parameter V > such that the quasi-interpolant with nodes in a neighborhood of [P, Q] 

Ml, v m = V-^ £ f(hm)r, (^=|P) 

d(hm,[P,Q])<r hVV 

approximates / with 

|/(x) - Ml v f{x)\ = 0((Vvh) N + e)\\f\\ w „ (1.4) 



y — hm\ 

dy 



for all x G [P,Q]. 
Then the integral 

£ M /(x) = MAW)(x) = ^" n/2 E / «a(x - y) r? f" 

d (fc m,[P,Q])<rW2> ^ V ^ 

gives a cubature of (1.1). 

Since KL\ is a bounded mapping between suitable function spaces, the differences /C^ ^/(x) — 
IC\f(x) behave like estimate (1.4). Therefore, to construct high order cubature formulas for 
(1.1), it remains to compute the integrals 

/ k a( - y= -y)v(y)dy 

J[p,Q] hVV 



for nodes with d(hm,[P,Q]) < rh\fT>. This is performed by using one-dimensional integral 
representations. As basis functions we take the tensor products of univariate basis functions 

2M 

5y(x) = HmMixj); m M {x 3 ) = ^ 2 L^ 2 \{x])^ (1.5) 

which satisfies the moment condition (1.2) of order N = 2M (cf. [11]), where are the 
generalized Laguerre polynomials 



4^^)Vv-). 7>-L 
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Using the representation with a tensor product integrand 

/ /C A (x - y)e" |y|2 dy = 
J[P,Q] 

\ /.-*"n^V - ^)) - erfc (f±±( Qj - i)))* (1,, 



we derive a tensor product representation of the integral operator which admits efficient cubature 
procedures for densities with separated approximation (Section 2). We will consider quasi- 
interpolants (2.1) on anisotropic grids which use different step size hj > 0,j = l,...,n along 
different space dimensions. If hj = rh , < r < 1 the error of the quasi-interpolant (2.1) is 
always 0{h N ). In Section 3 we provide numerical tests, showing that these formulas are accurate 
and provide approximation of order 0(h 6 ) up to dimension 10 s . 



2 Higher order cubature formula based on (1.6) 

In this section we describe a high order cubature of K\f in the case of rectangular domain in 
R n . Let 

n 

[P, Q] = {x = (xi, . . . , x n ) : Pj < Xj < Qj,j = 1, n} = JJ[P,-, Qj] . 

As basis functions we use (1.5). 

In order to apply also quasi-interpolants on rectangular grids {h\m\, . . . ,h n m n ),hj > 0, 
shortly denoted by {hm}, 



MwfW = V-^ Y, /( hm ) II *™ [ X \ h S 3 ) > (2-1) 

m£Z™ j=l ' 7 



we define the basis function ry(x) = YlV2M(ajXj), a,j > 0, and look for integral representations 
of the solution of 

n 

(-A + A 2 ) u = Yl X( Pj ,q 3 )(xj) rj 2 M(ajXj) . (2.2) 

3=1 

Here X(p j ,q j ) is t ne characteristic function of the interval (j)j,qj) with — oo < pj < qj < +oo, 
j = 1, . . . ,n. 

Theorem 2.1. Let Re A 2 > and n > 3. The solution of equation (2.2) in W l can be expressed 
by the one- dimensional integral 

I 7 2 

u(x) = j e~ A2 * /4 (^ M (ajXj,a 2 jt,ajPj) - ^ M (ajXj, a 2 jt, ajqjjj dt (2.3) 
o i= l 

where the function <3?m is given by 

e ~x 2 /(l+t) . e -F 2 (t,x,p) . 

$ M {x,t,p) = — -^j= — (erfc (F(t,x,p))V M (t,x) -j= QM{t,x,p)j 
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with the function 



F(t, X ,y) = ] J—(y- T ^), 
and Vm, Qm are polynomials in x of degree 2M — 2 and 2M — 3, respectively: 



Af-1 



^)-£(TT^- i/2, (iy- 

Q M (t, x, y) = 2 £ y- £ ^(y)^-i 



fc=l €=1 

2k\ fx \H e -i(F(t,x,y)) 



VT+tJ (i + t)fe+i/2 



//Re A 2 > 0, i/ien the representation (2.3) is valid for all n > 1. 
i??/ i/fc we denote the Hermite polynomials 

H k {x) = {-lfe* 2 ^- k e-* 2 . (2.4) 

Proof. The solution of (2.2) can be obtained explicitly by using the parabolic equation 

a(W-Aio + A 2 i« = 0, i > , (2.5) 

with the initial condition 

n 

w(x, 0) = Yl X(p j7qj )(xj) rj 2 M(a jXj ) . 
i=i 

Integrating (2.5) in i we derive 

T 

w(x,T) — io(x, 0) - (A - A 2 ) y w{x,t)dt = 0, 

o 

hence the solution of (2.2) is expressed as the one-dimensional integral 

oo 

u(x) = J w(x, t) dt 


provided it exists. Obviously, if w solves (2.5), then z = w e a2 * is the solution of the initial value 
problem for the heat equation 

n 

d t z - Az = , z(x, 0) = Yl Xip^q^ixj) V2M(ajXj) , 
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which has, by Poisson's formula, the solution 

. n 



(4?rt) 



1 f e -(oj*j-w) 2 /(4«i*) fhwru-i 

1 7rV2(4aft)V2 7 6 V2M{V 3 )dy 3 
where II (Pj;??) is the Cartesian product of the intervals (j>j,qj). Denoting 

oo 

<MM,p) = 4= / z- {x - y?/t mM{y)dy 

p 

we get the one-dimensional integral representation (2.3) of the solution of (2.2), provided this 
integral exists. Denoting 



oo 

j2k 



/d 
e ~ {x ~ y)2/t dj/ik e ~ y d y 



P 

and using the general representation [11, p. 55] 

r ?2M (x) = ^/ 2 ^ 1 ^A^e-N 2 , 

3=0 h 

we have 

k=0 



From 



Mx,t,P) = je-^^e-y' 2 dy = ^^L_ e ~ x2 ^ erfc (F(t,x,p)) , 



p 

for k > 1, integration by parts leads to 

2k-l 



w(x , , P ) = 5!!^ _ £ (-1)^ e-<-)V e- 



^=0 

and the definition (2.4) gives 



y=p 



d 2k- 



-v 



1 



dy 2k-£-l 

^ e -MV« _ (-!)< e-(*-vr* (y-x 
dy^ ~ UP £ \ Vt 
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In view of 



one gets for i < 2k 



^_erfc(x) = -^(-l/e- ;E2 ^_ 1 ( a ;), £>1, 



Q2k-£ r_]\2k-l r (fk-£ 

^ s - z erfc(F(t,x,p)) = (t(1 + t))fc _ </2 [d^ erfc(; 



z=F(t,x,p) 



n p -F 2 (t,x,p) 



H 2 k-e-i{F(t,x,p)) 



Therefore, since 



we obtain 



^^(M,*) = — (i + t)fc+1/2 g 2fc ^-^=j erfc(F(t,x,p)) 
V* e-^/Ci+t) e-P a (*.^) ^ ( 2k\ (-1) 



(1 + t)*+V2 

Thus simple transformations give 

¥>*(*,«, J>) =e- 2 /(l«) (erfc(F(t,i,p))/f 2t ( 



X 



'7rt 



VT+7/ 2(1 + t)*+V2 



^-l)/2 



= 1 



( 2 *W-<(: 



x \H^ l {F{t,x,p)) 



-Hp 



■ VT+t) {l + t) k + 1 / 2 
Using the relation H 2k (x) = (-l) fe 4 fe fe! L<T 1/2) ( X 2 ) we find therefore 



p — x 



H 2 k-i 



e -x 2 /(i+t) erf c (F(i,a;,p)) 



M-l 1 o 

V 1 T {-m( x \ 

(l +t )fc+l/2 k \l + t ) 



fc=0 ' €=1 



7T 



e -x 2 /(l+t) e -F 2 (t,x,p) 

y=— y erfc (F(t, x,p))V M (t,x) -= Qm {t,x,p)j. 



2^ 



The polynomials Vm^iX) and Q,M(t,x,p) for M = 1,2,3 are given by 
Pi(t,x) = - , P 2 (t,x) = Pi(t,x) + 



(1 + t)V2 ' ' ' ' iV ' ' 2(l+t)3/ 2 (l + t) 5 /2' 

. . . . 3 3 x x 

V 3 (t, x) = V 2 (t, x) + — — - - — + 



8(l + t) 5 / 2 2(l + t) 7 /2 2(l+i) 9 /2 ' 

Qi(t,x,p) = 0, Q 2 (t,x,p) = ^^(^ 

/ 2x 3 2px 2 - 5x (2p 2 - 5)x - 3p , 



Remark 2.1. Smce /or positive r 

< erfc(r) < e -7 * 2 and 2 - e" r2 < erfc(-r) < 2 



/rom i/ie relation 



F (t,x,p) =p + — 



t l + t 

we get 

| e -*7(i+*) erfc (F(t, x, p)) | < e~ p2 »/ p > 

and 

| e -* 2 /(i+t) erfc (F(t, x, p)) - 2 e - a,a /( 1+ *) | < e~ p2 if p<0. 
Thus for sufficiently large \p\ 

7r-V2 e -*V(i+t)p M (i,z ) + 0(e-P 2 ) z /p<0, 
$M(X ''' P) -\ 0(e"f 2 ) z/p>0, 

and therefore, for sufficiently large r one can use the approximation 

* I + \ * < * \ ~ / °' p,q>rorp,q<-r, 

with the error 0(e~ r ). Similarly, if q — p > 2r, then 

it, ( + \ rt, ( + \ ~ / ®M(x,t,p), -r<p<r, 

*M(x,t, P )-*M(x,t, q )~[ n - 1/2 e -W + t) pM{tiX) _t M{Xitiq)i _ r<0<r . 

3 Implementation and numerical results 

We compute the cubature formula 



/c A ,h/(x) = v~ n i 2 Y, /( hm ) / Kx(x-y)f[m M ( 



hjVV 
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where f2 r h = 11^=1 (-^j ~~ rhj^/V, Qj + rhj^/V), using the tensor product representation of The- 
orem 2.1. At the grid points hk = (h\ki, . . . , h n k n ) we obtain 



oo 

2 



e -A»t/4 



[P,Q] 1=1 



n( / Jx, j Tfl j t Pj h j Tfl j / k j Tfl j t j h j Tfl j \ \ 

and therefore 

/C A ,h/(x)= Yl /( hm ) b S> 

hm£n r h 

where we introduce the one-dimensional integral 

00 n 

bg£ = ^ / e" A2</4 II fa. ^) " 6 U Wi)) * ( 3 - 2 ) 
o -? =1 

and use the abbreviation 

b j (P) - e -(k- m f/(v(i + t)) ( eifc(F( _L_ p ~ hjm ))v M ( f fc ~ m ) 

/ , t k — m P — Lm,\ , t k — m P — hjm\\ ,, ^ 



-7r x / 2 exp 



According to Remark 2.1, for appropriately chosen r > we can set within a given accuracy 
K,m(P) = <-m = -~ 1/2 e-^) 2 ^(i +t )) V M (^ ^ if P - h 3 m < -rh.Vv, 

K,m( P ) =° if P " ^ m > rh 3 > 

which speeds up the computation of (3.2). In particular, we can split (3.1) into 

4 M h V(hk)= E/( hm )4_L+ £ 7(hm)bg2. (3.3) 

hmeC rh hmeSJ rh \Jl rh 

where f2 r h = Ilj=i(-fj + rhjVV, Qj — rhjVV), and the coefficients in the convolutional sum are 
given by 



1 — -7 j 1 • I T a?, dt 

k 4X)™/ 2 y 11 fc j 

J =1 

00 
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Following [12] the one-dimensional integrals of and are transformed to integrals 

over R with integrands decaying doubly exponentially by making the substitutions 

t = e«, £ = a(<r + e CT ), a = p[u - e" M ) (3.4) 

with certain positive constants a, (3, and the computation is based on the classical trapezoidal 
rule. Then the tensor product structure of the integrands allows the efficient computation of 
the coefficients bj^ and aj^. Moreover, the computation of the convolutional sum is very 
efficient for integrands, which allow a separated representation, i.e., for given accuracy e they 
can be represented as a sum of products of vectors in dimension 1 

R n 

/(/iimi, . . . , h m m n ) = Y[ fj P \ h j m j) + C( e ) • 

p=i j=i 

In [7] we have described this approach to the fast computation of high dimensional volume 
potentials for compactly supported integrands. To compute the convolutional sum 

E 4 M L/(hm) 

hm£f2 r h 

we get after the substitutions 

j—i j 

where we set 



<&(u) = exp(a/3(ti — exp(— u)) + a exp(/3(u — exp(— u)))), 
= &(u)ap(l + e~ M )(l + exp(/3(u - exp(-u)))). 

The quadrature with the trapezoidal rule with step size r 



( M ) ~ T V C - X2 ^ s c -!k| 2 /(P(l+$(s r))) TT j, M ST ) Ji_\ 



s=-AT i=l J 

provides the approximation via one-dimensional discrete convolutions 



hmeSi rh V ' p=l s=-N 



We provide some numerical tests to the approximation of the potential /C\/ over the cube 
[— 1, l] n , n > 3, with the density 



n n 

tip), 



/(x) = (-a + a 2 ) n «(^-) = e n (^). x = - . *») e ^ ( 3 - s ) 

jr' = l p=l j=l 
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&\x) = u(x) if j ± p; ff\x) = -u"(x) + —u(x) if j = p. 



n 



Let be an extension of outside the interval [—1, 1] with preserved smoothness and 



p=i j=i 



By using Hestenes reflection principle ( [5] ) we construct an extension of outside the interval 
[-1,1] as 

AT+l 

£c,ff\-a,(x + l)-l), x<-l 

s=l 



-1< X < 1 



N+l 



^c s /j p) (-a s (x-l) + l), x>l 



s=l 



where ai,...,ajv+i are different positive constants and the coefficients c^r = {c±, cat+i} are 
the unique solution of the (N + 1) x (N + 1) system of linear equations 



N+l 
s=l 



1, k = 0,...,N. 



~<n) (p) 

We provide results for /j = /j and three different Hestenes extensions corresponding to 
a s = 2~ s (Extension 1) , a s = s" 1 (Extension 2), a s = s (Extension 3). 

The approximation values are computed by the cubature formula (3.3) for hj = h, j = 1, n. 
To have the saturation error comparable with the double precision rounding errors, we have 
chosen the parameter T> = 4. 

In Tables 1, 2 and 3 we report on the absolute error and the approximation rate for the 
three-dimensional potential JC\f, when u(x) = cos 2 (7rx/2) (Table 1), u(x) = (x 2 — l) 3 (Table 2) 
and u(x) = (x 2 — l) 2 (Table 3), in the case A 2 = 1 and A 2 = We have chosen the parameters 
a = 2, (3 = 2 in the transformations (3.4) and r = 0.005, N\ = —Nq = 300 in the quadrature 
formula. The numerical results confirm the h 2 — , h 4 — and, respectively, h 6 — convergence of the 
cubature formulas (3.3) when M = 1,2,3. For extensions 1, 2 and 3 the numerical results are 

~fp) (n) 

similar with those if using 

H = H ■ In Table 3 

we see that the error of the approximate 
quasi-interpolant of order 6 has reached the saturation bound. This is a feature of the method 
that approximate quasi-interpolant of order N reproduces polynomials of degree < N up to the 
saturation error. 

To check the effectiveness of the method for very high dimension n we computed the potential 
over [-1, l] n of the density (3.5) with u(x) = 1 - sin(7rx 2 /2) (Table 4) and u(x) = e x (l - x 2 ) 2 



(Table 5) in dimension n = 10*, i 
b = 5, r = 0.003, N = -40, iVi = 
predicted approximation rate 6, also for very large n and the error scales linearly in the space 
dimension. 



1,..,8 and different extensions. We have chosen a = 6, 
200. The results show that JC^\ approximates with the 
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A 2 = l: 







M = 


1 


M = 


2 


M = 


3 


/(*) 


h- 1 


error 


rate 


error 


rate 


error 


rate 




10 


0.822E-01 




0.414E-02 




0.135E-03 






20 


0.219E-01 


1.9062 


0.272E-03 


3.9267 


0.223E-05 


5.9201 


/(*) 


40 


0.557E-02 


1.9760 


0.172E-04 


3.9821 


0.354E-07 


5.9800 




80 


0.140E-02 


1.9940 


0.108E-05 


3.9955 


0.555E-09 


5.9950 




160 


0.350E-03 


1.9985 


0.675E-07 


3.9989 


0.867E-11 


5.9987 




320 


0.875E-04 


1.9996 


0.422E-08 


3.9997 


0.136E-12 


5.9994 




10 


0.821E-01 




0.413E-02 




0.135E-03 






20 


0.219E-01 


1.9057 


0.272E-03 


3.9265 


0.223E-05 


5.9201 


cxt 1 


40 


0.557E-02 


1.9760 


0.172E-04 


3.9820 


0.354E-07 


5.9800 




80 


0.140E-02 


1.9940 


0.108E-05 


3.9955 


0.554E-09 


5.9961 




160 


0.350E-03 


1.9985 


0.675E-07 


3.9989 


0.825E-11 


6.0692 




320 


0.875E-04 


1.9996 


0.422E-08 


3.9997 


0.789E-12 


3.3868 




10 


0.826E-01 




0.422E-02 




0.140E-03 






20 


0.219E-01 


1.9138 


0.273E-03 


3.9520 


0.224E-05 


5.9686 


cxt 2 


40 


0.557E-02 


1.9769 


0.172E-04 


3.9850 


0.354E-07 


5.9856 




80 


0.140E-02 


1.9941 


0.108E-05 


3.9959 


0.554E-09 


5.9967 




160 


0.350E-03 


1.9985 


0.675E-07 


3.9989 


0.883E-11 


5.9718 




320 


0.875E-04 


1.9996 


0.422E-08 


3.9997 


0.120E-12 


6.1971 




10 


0.946E-01 




0.139E-01 




0.260E-01 






20 


0.224E-01 


2.0769 


0.771E-03 


4.1771 


0.871E-04 


8.2194 


cxt 3 


40 


0.559E-02 


2.0047 


0.228E-04 


5.0788 


0.111E-05 


6.2957 




80 


0.140E-02 


1.9977 


0.113E-05 


4.3396 


0.341E-08 


8.3438 




160 


0.350E-03 


1.9990 


0.679E-07 


4.0529 


0.147E-10 


7.8633 




320 


0.875E-04 


1.9997 


0.422E-08 


4.0067 


0.147E-12 


6.6382 



A 2 = 1 + 







M = 


1 


M = 


2 


M = 


3 


/(x) 


h- 1 


error 


rate 


error 


rate 


error 


rate 




10 


0.815E-01 




0.410E-02 




0.134E-03 






20 


0.217E-01 


1.9060 


0.270E-03 


3.9267 


0.221E-05 


5.9201 


/(x) 


40 


0.553E-02 


1.9760 


0.171E-04 


3.9821 


0.351E-07 


5.9800 




80 


0.139E-02 


1.9940 


0.107E-05 


3.9955 


0.550E-09 


5.9950 




160 


0.347E-03 


1.9985 


0.669E-07 


3.9989 


0.860E-11 


5.9987 




320 


0.868E-04 


1.9996 


0.418E-08 


3.9997 


0.135E-12 


5.9974 




10 


0.814E-01 




0.410E-02 




0.134E-03 






20 


0.217E-01 


1.9055 


0.270E-03 


3.9265 


0.221E-05 


5.9201 


cxt 1 


40 


0.553E-02 


1.9759 


0.171E-04 


3.9820 


0.351E-07 


5.9800 




80 


0.139E-02 


1.9940 


0.107E-05 


3.9955 


0.550E-09 


5.9959 




160 


0.347E-03 


1.9985 


0.669E-07 


3.9989 


0.826E-11 


6.0555 




320 


0.868E-04 


1.9996 


0.418E-08 


3.9997 


0.710E-12 


3.5419 




10 


0.819E-01 




0.417E-02 




0.139E-03 






20 


0.218E-01 


1.9127 


0.270E-03 


3.9490 


0.222E-05 


5.9631 


cxt 2 


40 


0.553E-02 


1.9768 


0.171E-04 


3.9846 


0.351E-07 


5.9849 




80 


0.139E-02 


1.9941 


0.107E-05 


3.9959 


0.550E-09 


5.9964 




160 


0.347E-03 


1.9985 


0.669E-07 


3.9989 


0.873E-11 


5.9767 




320 


0.868E-04 


1.9996 


0.418E-08 


3.9997 


0.122E-12 


6.1594 




10 


0.924E-01 




0.130E-01 




0.238E-01 






20 


0.222E-01 


2.0586 


0.717E-03 


4.1823 


0.799E-04 


8.2188 


cxt 3 


40 


0.554E-02 


2.0011 


0.220E-04 


5.0283 


0.101E-05 


6.3010 




80 


0.139E-02 


1.9973 


0.111E-05 


4.3051 


0.313E-08 


8.3370 




160 


0.347E-03 


1.9989 


0.673E-07 


4.0461 


0.139E-10 


7.8181 




320 


0.868E-04 


1.9997 


0.419E-08 


4.0058 


0.145E-12 


6.5840 



Table 1: Absolute errors and approximation rates for /C\/(0.3, 0.3, 0) using JC\ . /(0.3, 0.3, 0) 
with the density / given in (3.5) with u(x) = cos 2 (ttx/2) and different extensions, M = 1,2,3, 
A 2 = 1 and A 2 = 1 + i. 
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A 2 = l: 







M = 


1 


M = 


2 


M = 


3 


/(*) 


h- 1 


error 


rate 


error 


rate 


error 


rate 




10 


0.673E-01 




0.626E-02 




0.427E-04 






20 


0.159E-01 


2.0819 


0.392E-03 


3.9965 


0.668E-06 


5.9997 


/(*) 


40 


0.391E-02 


2.0238 


0.246E-04 


3.9970 


0.104E-07 


6.0000 




80 


0.973E-03 


2.0062 


0.154E-05 


3.9991 


0.163E-09 


6.0000 




160 


0.243E-03 


2.0016 


0.960E-07 


3.9998 


0.255E-11 


6.0000 




320 


0.607E-04 


2.0004 


0.600E-08 


3.9999 


0.398E-13 


6.0002 




10 


0.637E-01 




0.634E-02 




0.427E-04 






20 


0.157E-01 


2.0254 


0.393E-03 


4.0094 


0.668E-06 


5.9997 


cxt 1 


40 


0.389E-02 


2.0075 


0.246E-04 


4.0003 


0.104E-07 


5.9995 




80 


0.972E-03 


2.0020 


0.154E-05 


3.9999 


0.156E-09 


6.0635 




160 


0.243E-03 


2.0005 


0.961E-07 


4.0000 


0.389E-11 


5.3255 




320 


0.607E-04 


2.0001 


0.600E-08 


4.0000 


0.603E-12 


2.6899 




10 


0.603E-01 




0.644E-02 




0.427E-04 






20 


0.154E-01 


1.9662 


0.395E-03 


4.0264 


0.668E-06 


5.9997 


cxt 2 


40 


0.388E-02 


1.9925 


0.246E-04 


4.0052 


0.104E-07 


6.0003 




80 


0.971E-03 


1.9983 


0.154E-05 


4.0012 


0.163E-09 


6.0019 




160 


0.243E-03 


1.9996 


0.961E-07 


4.0003 


0.224E-11 


6.1838 




320 


0.607E-04 


1.9999 


0.600E-08 


4.0001 


0.408E-12 


2.4557 




10 


0.291E-01 




0.626E-02 




0.427E-04 






20 


0.133E-01 


1.1335 


0.392E-03 


3.9965 


0.668E-06 


5.9997 


cxt 3 


40 


0.374E-02 


1.8264 


0.246E-04 


3.9970 


0.104E-07 


6.0000 




80 


0.963E-03 


1.9586 


0.154E-05 


3.9991 


0.163E-09 


6.0000 




160 


0.224E-03 


1.9894 


0.960E-07 


3.9998 


0.255E-11 


6.0000 




320 


0.607E-04 


1.9975 


0.600E-08 


3.9999 


0.398E-13 


6.0001 



A 2 = 1 + 







M = 


1 


M = 


2 


M = 


3 


/(x) 


h- 1 


error 


rate 


error 


rate 


error 


rate 




10 


0.604E-01 




0.572E-02 




0.441E-04 






20 


0.142E-01 


2.0834 


0.358E-03 


3.9963 


0.690E-06 


5.9997 


/(x) 


40 


0.350E-02 


2.0242 


0.224E-04 


3.9969 


0.108E-07 


6.0000 




80 


0.872E-03 


2.0062 


0.140E-05 


3.9991 


0.168E-09 


6.0000 




160 


0.218E-03 


2.0016 


0.878E-07 


3.9998 


0.263E-11 


6.0000 




320 


0.544E-04 


2.0004 


0.548E-08 


3.9999 


0.410E-13 


6.0025 




10 


0.572E-01 




0.579E-02 




0.441E-04 






20 


0.140E-01 


2.0271 


0.360E-03 


4.0096 


0.690E-06 


5.9997 


cxt 1 


40 


0.349E-02 


2.0080 


0.225E-04 


4.0004 


0.108E-07 


5.9996 




80 


0.871E-03 


2.0021 


0.140E-05 


4.0000 


0.163E-09 


6.0465 




160 


0.218E-03 


2.0006 


0.878E-07 


4.0000 


0.372E-11 


5.4561 




320 


0.544E-04 


2.0001 


0.548E-08 


4.0000 


0.539E-12 


2.7853 




10 


0.542E-01 




0.589E-02 




0.441E-04 






20 


0.138E-01 


1.9681 


0.361E-03 


4.0272 


0.690E-06 


5.9997 


cxt 2 


40 


0.348E-02 


1.9931 


0.225E-04 


4.0055 


0.108E-07 


6.0002 




80 


0.870E-03 


1.9984 


0.140E-05 


4.0013 


0.168E-09 


6.0014 




160 


0.218E-03 


1.9996 


0.878E-07 


4.0003 


0.240E-11 


6.1310 




320 


0.544E-04 


1.9999 


0.548E-08 


4.0001 


0.365E-12 


2.7174 




10 


0.261E-01 




0.803E-02 




0.441E-04 






20 


0.119E-01 


1.1338 


0.560E-03 


3.8421 


0.690E-06 


5.9997 


cxt 3 


40 


0.335E-02 


1.8275 


0.268E-04 


4.3875 


0.108E-07 


6.0000 




80 


0.862E-03 


1.9590 


0.148E-05 


4.1767 


0.168E-09 


6.0000 




160 


0.217E-03 


1.9899 


0.890E-07 


4.0553 


0.263E-11 


6.0000 




320 


0.544E-04 


1.9975 


0.550E-08 


4.0151 


0.410E-13 


6.0030 



Table 2: Absolute errors and approximation rates for /C\/(0.5, 0.5, 0.5) using /C^ h /(0.5, 0.5, 0.5) 
with the density / given in (3.5) with u(x) = (x 2 — l) 3 and different extensions, M = 1,2,3, 
A 2 = 1 and A 2 = 1 + i. 
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M = 


1 


M = 


2 


M = 


3 


/(x) 


h- 1 


error 


rate 


error 


rate 


error 


rate 




10 


0.935E-01 




0.166E-02 




0.222E-15 






20 


0.241E-01 


1.9564 


0.104E-03 


3.9984 


0.777E-15 




/(*) 


40 


0.607E-02 


1.9883 


0.647E-05 


3.9999 


0.111E-15 






80 


0.152E-02 


1.9970 


0.405E-06 


4.0000 


0.555E-16 






160 


0.380E-03 


1.9993 


0.253E-07 


4.0000 


0.555E-16 






320 


0.951E-04 


1.9998 


0.158E-08 


4.0000 


0.222E-15 






10 


0.941E-01 




0.166E-02 




0.779E-10 






20 


0.241E-01 


1.9632 


0.104E-03 


3.9984 


0.336E-10 


1.2133 


cxt 1 


40 


0.607E-02 


1.9903 


0.647E-05 


3.9999 


0.143E-10 


1.2318 




80 


0.152E-02 


1.9975 


0.405E-06 


4.0000 


0.628E-11 


1.1873 




160 


0.380E-03 


1.9994 


0.253E-07 


4.0000 


0.160E-12 


5.2966 




320 


0.951E-04 


1.9998 


0.158E-08 


4.0000 


0.268E-12 


-0.7471 




10 


0.946E-01 




0.166E-02 




0.201E-11 






20 


0.242E-01 


1.9684 


0.104E-03 


3.9984 


0.133E-11 


0.5895 


ext 2 


40 


0.607E-02 


1.9920 


0.647E-05 


3.9999 


0.139E-11 


-0.0557 




80 


0.152E-02 


1.9980 


0.405E-06 


4.0000 


0.641E-13 


4.4336 




160 


0.380E-03 


1.9995 


0.253E-07 


4.0000 


0.532E-12 


-3.0524 




320 


0.951E-04 


1.9999 


0.158E-08 


4.0000 


0.404E-12 


0.3980 




10 


0.983E-01 




0.166E-02 




0.222E-15 






20 


0.245E-01 


2.0041 


0.104E-03 


3.9984 


0.722E-15 




ext 3 


40 


0.610E-02 


2.0066 


0.647E-05 


3.9999 


0.111E-15 






80 


0.152E-02 


2.0022 


0.405E-06 


4.0000 


0.555E-16 






160 


0.380E-03 


2.0006 


0.253E-07 


4.0000 


0.555E-16 






320 


0.951E-04 


2.0002 


0.158E-08 


4.0000 


0.111E-15 





A 2 = 1 + 







M = 


1 


M = 


2 


M = 


3 


/(x) 


h- 1 


error 


rate 


error 


rate 


error 


rate 




10 


0.869E-01 




0.168E-02 




0.220E-14 






20 


0.224E-01 


1.9541 


0.105E-03 


3.9983 


0.729E-15 




/(x) 


40 


0.565E-02 


1.9878 


0.655E-05 


3.9999 


0.397E-15 




80 


0.142E-02 


1.9969 


0.410E-06 


4.0000 


0.555E-16 






160 


0.354E-03 


1.9992 


0.256E-07 


4.0000 


0.128E-15 






320 


0.886E-04 


1.9998 


0.160E-08 


4.0000 


0.906E-16 






10 


0.875E-01 




0.168E-02 




0.695E-10 






20 


0.225E-01 


1.9617 


0.105E-03 


3.9983 


0.301E-10 


1.2058 


cxt 1 


40 


0.566E-02 


1.9899 


0.655E-05 


3.9999 


0.128E-10 


1.2359 




80 


0.142E-02 


1.9974 


0.410E-06 


4.0000 


0.563E-11 


1.1851 




160 


0.354E-03 


1.9994 


0.256E-07 


4.0000 


0.144E-12 


5.2918 




320 


0.886E-04 


1.9998 


0.160E-08 


4.0000 


0.240E-12 


-0.7425 




10 


0.880E-01 




0.168E-02 




0.179E-11 






20 


0.225E-01 


1.9675 


0.105E-03 


3.9983 


0.119E-11 


0.5930 


cxt 2 


40 


0.566E-02 


1.9917 


0.655E-05 


3.9999 


0.124E-11 


-0.0642 




80 


0.142E-02 


1.9979 


0.410E-06 


4.0000 


0.577E-13 


4.4276 




160 


0.354E-03 


1.9995 


0.256E-07 


4.0000 


0.476E-12 


-3.0451 




320 


0.886E-04 


1.9999 


0.160E-08 


4.0000 


0.361E-12 


0.3971 




10 


0.918E-01 




0.168E-02 




0.215E-14 






20 


0.228E-01 


2.0074 


0.105E-03 


3.9983 


0.625E-15 




ext 3 


40 


0.568E-02 


2.0074 


0.655E-05 


3.9999 


0.296E-15 






80 


0.142E-02 


2.0024 


0.410E-06 


4.0000 


0.706E-17 






160 


0.354E-03 


2.0006 


0.256E-07 


4.0000 


0.794E-16 






320 


0.886E-04 


2.0002 


0.160E-08 


4.0000 


0.119E-15 





Table 3: Absolute errors and approximation rates for JC\f(0A, 0.5, 0) using JC\ . /(0.4, 0.5, 0) 
with the density / given in (3.5) with u(x) = (1 — x 2 ) 2 and different extensions, with M = 1,2, 3, 
A 2 = 1 and A 2 = 1 + i. 
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/(x) 


?i 


10 


10 2 


10 3 


10 4 




ft- 1 


error rate 


error rate 


error rate 


error rate 


/(*) 


10 

20 
40 
80 
160 
320 


0.338E-03 
0.605E-05 5.8020 
0.976E-07 5.9541 
0.154E-08 5.9887 
0.241E-10 5.9971 
0.376E-12 5.9982 


0.459E-02 
0.732E-04 5.9727 
0.115E-05 5.9966 
0.179E-07 5.9994 
0.280E-09 6.0013 
0.513E-11 5.7677 


0.487E-01 

0.746E-03 6.0282 
0.117E-04 5.9991 
0.182E-06 5.9999 
0.285E-08 6.0000 
0.445E-10 6.0005 


0.703E+00 
0.751E-02 6.5491 
0.117E-03 6.0070 
0.183E-05 6.0000 
0.285E-07 5.9999 
0.446E-09 5.9985 


/(x) 


n 


10 5 


10 6 


10 7 


10 8 






error rate 


error rate 


error rate 


error rate 


/(x) 


20 
40 
80 
160 
320 


0.794E-01 
0.117E-02 6.0852 
0.183E-04 6.0012 
0.285E-06 5.9992 
0.451E-08 5.9842 


0.145E+01 
0.118E-01 6.9443 
0.183E-03 6.0133 
0.286E-05 5.9975 
0.478E-07 5.9030 


0.129E+00 
0.183E-02 6.1364 
0.286E-04 5.9985 
0.510E-06 5.8096 


0.348E+01 
0.185E-01 7.5527 
0.286E-03 6.0174 
0.517E-05 5.7889 




/(x) 


n 


10 


10 2 


10 3 


10 4 






error rate 


error rate 


error rate 


error rate 


ext 
1 


10 

20 
40 
80 
160 
320 


0.352E-03 
0.611E-05 5.8476 
0.978E-07 5.9652 
0.154E-08 5.9892 
0.230E-10 6.0635 
0.650E-12 5.1472 


0.459E-02 
0.732E-04 5.9726 
0.115E-05 5.9966 
0.179E-07 5.9994 
0.280E-09 6.0013 
0.513E-11 5.7677 


0.487E-01 

0.746E-03 6.0282 
0.117E-04 5.9991 
0.182E-06 5.9999 
0.285E-08 6.0000 
0.445E-10 6.0005 


0.703E+00 
0.751E-02 6.5491 
0.117E-03 6.0070 
0.183E-05 6.0000 
0.285E-07 5.9999 
0.446E-09 5.9985 


/(x) 


n 


10 5 


10 (i 


10 7 


10 8 






error rate 


error rate 


error rate 


error rate 


ext 
1 


20 
40 
80 
160 
320 


0.794E-01 
0.117E-02 6.0852 
0.183E-04 6.0012 
0.285E-06 5.9992 
0.451E-08 5.9842 


0.145E+01 
0.118E-01 6.9443 
0.183E-03 6.0133 
0.286E-05 5.9975 
0.478E-07 5.9030 


0.129E+00 
0.183E-02 6.1364 
0.286E-04 5.9985 
0.510E-06 5.8096 


0.348E+01 
0.185E-01 7.5527 
0.286E-03 6.0174 
0.517E-05 5.7889 




/(x) 


n 


10 


10 2 


10 3 


10 4 






error rate 


error rate 


error rate 


error rate 


ext 
2 


10 

20 
40 
80 
160 
320 


0.415E-03 
0.632E-05 6.0374 
0.985E-07 6.0037 
0.154E-08 5.9994 
0.241E-10 5.9999 
0.408E-12 5.8832 


0.459E-02 
0.732E-04 5.9727 
0.115E-05 5.9966 
0.179E-07 5.9994 
0.280E-09 6.0013 
0.513E-11 5.7677 


0.487E-01 

0.746E-03 6.0282 
0.117E-04 5.9991 
0.182E-06 5.9999 
0.285E-08 6.0000 
0.445E-10 6.0005 


0.703E+00 
0.751E-02 6.5491 
0.117E-03 6.0070 
0.183E-05 6.0000 
0.285E-07 5.9999 
0.446E-09 5.9985 


/(x) 


n 


10 5 


10 6 


10 7 


10 8 




h~ L 


error rate 


error rate 


error rate 


error rate 


ext 
2 


20 
40 
80 
160 
320 


0.794E-01 
0.117E-02 6.0852 
0.183E-04 6.0012 
0.285E-06 5.9992 
0.451E-08 5.9842 


0.145E+01 
0.118E-01 6.9443 
0.183E-03 6.0133 
0.286E-05 5.9975 
0.478E-07 5.9030 


0.129E+00 
0.183E-02 6.1364 
0.286E-04 5.9985 
0.510E-06 5.8096 


0.348E+01 
0.185E-01 7.5527 
0.286E-03 6.0174 
0.517E-05 5.7889 



Table 4: Absolute errors and approximation rates for /Ca/(0.5, 0..., 0) using JC X h f (0.5,0, . 
with the density / given in (3.5) with u{x) = 1 — sin(-7ra; 2 /2) and different extensions , n = 
i = 1,...,8, A 2 = 1. 
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/(x) 


n 


10 


10 2 


10 3 


10 4 




ft- 1 


error rate 


error rate 


error rate 


error rate 


/(x) 


10 

20 
40 
80 
160 
320 


0.699E-03 
0.106E-04 6.0400 
0.165E-06 6.0100 
0.257E-08 6.0025 
0.402E-10 6.0005 
0.632E-12 5.9909 


0.596E-02 
0.902E-04 6.0453 
0.140E-05 6.0105 
0.218E-07 6.0026 
0.341E-09 6.0017 
0.491E-11 6.1156 


0.595E-01 

0.880E-03 6.0792 
0.136E-04 6.0111 
0.213E-06 6.0026 
0.332E-08 6.0006 
0.585E-10 5.9998 


0.759E+00 
0.881E-02 6.4288 
0.136E-03 6.0162 
0.212E-05 6.0027 
0.332E-07 6.0005 
0.519E-09 5.9973 


/(x) 


n 


10 5 


10 6 


10 7 


10 8 




h- L 


error rate 


error rate 


error rate 


error rate 


/(x) 


20 
40 
80 
160 
320 


0.913E-01 
0.136E-02 6.0671 
0.212E-04 6.0035 
0.332E-06 5.9994 
0.526E-08 5.9779 


0.134E+01 
0.137E-01 6.6101 
0.212E-03 6.0113 
0.332E-05 5.9966 
0.572E-07 5.8594 


0.145E+00 
0.212E-02 6.0906 
0.333E-04 5.9966 
0.632E-06 5.7186 


0.267E+01 
0.214E-01 6.9639 
0.333E-03 6.0087 
0.646E-05 5.6865 




/(x) 


n 


10 


10 2 


10 3 


10 4 






error rate 


error rate 


error rate 


error rate 


ext 
2 


10 

20 
40 
80 
160 
320 


0.690E-03 
0.106E-04 6.0254 
0.165E-06 6.0068 
0.257E-08 6.0019 
0.401E-10 6.0046 
0.676E-12 5.8884 


0.596E-02 
0.902E-04 6.0453 
0.140E-05 6.0105 
0.218E-07 6.0026 
0.341E-09 6.0017 
0.491E-11 6.1156 


0.595E-01 

0.880E-03 6.0792 
0.136E-04 6.0111 
0.213E-06 6.0026 
0.332E-08 6.0006 
0.519E-10 5.9998 


0.759E+00 
0.881E-02 6.4288 
0.136E-03 6.0162 
0.212E-05 6.0027 
0.332E-07 6.0005 
0.519E-09 5.9973 


/(x) 


n 


10 5 


10 6 


10 7 


10 8 




h- L 


error rate 


error rate 


error rate 


error rate 


ext 
2 


20 
40 
80 
160 
320 


0.913E-01 
0.136E-02 6.0671 
0.212E-04 6.0035 
0.332E-06 5.9994 
0.526E-08 5.9779 


0.134E+01 
0.137E-01 6.6101 
0.212E-03 6.0113 
0.332E-05 5.9966 
0.572E-07 5.8594 


0.145E+00 
0.212E-02 6.0906 
0.333E-04 5.9966 
0.632E-06 5.7186 


0.267E+01 
0.214E-01 6.9639 
0.333E-03 6.0087 
0.646E-05 5.6865 




/(x) 


n 


10 


10 2 


10 3 


10 4 






error rate 


error rate 


error rate 


error rate 


ext 
3 


10 

20 
40 
80 
160 
320 


0.156E-01 
0.165E-04 9.8811 
0.943E-07 7.4538 
0.110E-08 6.4188 
0.333E-10 5.0496 
0.602E-12 5.7901 


0.590E-02 
0.901E-04 6.0349 
0.140E-05 6.0091 
0.218E-07 6.0021 
0.340E-09 6.0016 
0.491E-11 6.1156 


0.595E-01 
0.880E-03 6.0791 
0.136E-04 6.0111 
0.213E-06 6.0026 
0.332E-08 6.0006 
0.519E-10 5.9998 


0.759E+00 
0.881E-02 6.4288 
0.136E-03 6.0162 
0.212E-05 6.0027 
0.332E-07 6.0005 
0.519E-09 5.9973 


/(x) 


n 


10 5 


10 6 


10 7 


10 8 






error rate 


error rate 


error rate 


error rate 


ext 
3 


20 
40 
80 
160 
320 


0.913E-01 
0.136E-02 6.0671 
0.212E-04 6.0035 
0.332E-06 5.9994 
0.526E-08 5.9779 


0.134E+01 
0.137E-01 6.6101 
0.212E-03 6.0113 
0.332E-05 5.9966 
0.572E-07 5.8599 


0.145E+00 
0.212E-02 6.0906 
0.333E-04 5.9966 
0.632E-06 5.7186 


0.267E+01 
0.214E-01 6.9639 
0.333E-03 6.0087 
0.646E-05 5.6865 



Table 5: Absolute errors and approximation rates for K.\f(0A, 0.4, 0, 0) using 
/C^/(0.4,0.4,0,...,0) with the density / given in (3.5) with u(x) = e x (l - x 2 ) 2 and differ- 
ent extensions , n = 10*, i = 1, 8, A 2 = 1. 



17 



